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Abstract 

We study the most general, relativistic, constituent qq meson wave function 
within a new covariant framework. We find that by including a tensor wave 
function component, a pure valence quark model is now capable of reproduc- 
ing not only all static pion data (/„-, (r%)) but also the distribution amplitude, 
form factor (F 7T (Q 2 )), and structure functions. Further, our generalized spin 
wave function provides a much better detailed description of meson prop- 
erties than models using a simple relativistic extension of the S = L = 

nonrelativistic wave function. 
12.40.Aa, 14.40.Aq, 13.60.-r 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Since an exact solution to a bound state problem in QCD is still unavailable many 
approximate treatments have been developed. Among them, the constituent quark model 
has perhaps received the most attention and is widely regarded as a very efficient and effective 
tool in description of hadronic phenomena (T|-§[. The complex low energy structure of QCD 
currently precludes an unambiguous identification of the complete degrees of freedom and it 
is therefore important to continue to advance, refine and test the valence quark dominance 
approximation. Prospects are encouraging that new, precision data provided by future 
CEBAF experiments will significantly clarify this situation and also detail the role of exotic 
quark and/or gluon configurations. However, before the relative importance of valence 
versus exotic configurations can be established, uncertainties in the description of hadronic 
amplitudes due to valence quark model approximations must be reduced. In a previous 
study |4| we investigated alternative model approaches by comparing different relativistic 
formulations for the light pseudoscalar mesons. In particular, we developed and examined 
a covariant variable front approach which permitted quantitatively assessing the relevance 
of Lorentz covariance for any constituent formulation. We also established that pion and 
the kaon static properties can be well described by relativistic models utilizing constituent 
qq meson wave functions represented by the product of a noninteracting spinor component 
and a momentum space orbital amplitude. It is believed that such models will be able to 
describe low energy mesonic data to within a 10-20%. Unfortunately, there is still a large 
energy gap between this region and the energy-momentum scale where asymptotic freedom 
dominates || and a clear need exists for an improved, QCD related quark model to describe 
this intermediated regime. This is especially true for analyzing electromagnetic form factors 
with Q 2 > 2 — 3 GeV 2 for which all quark models described in Ref. |4j] fail to describe 
the data. As mentioned above, however, before making any exotic extensions of the quark 
model an improved framework for the qq system must be developed. The purpose of the 
present paper is to report one such attempt which utilizes a more general quark meson wave 
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functions. 

In this article we extend our covariant variable front quark model by incorporating a 
more general spinor wave function with tensor components. Although not an observable, 
the wave function is constrained by results from QCD studies of moments of the distribution 
amplitudes H-H • Additional information comes from meson structure functions measured in 



Drell-Yan experiments fllOHl^l. In the next section we briefly review the basic assumptions 
of the covariant quark model, and detail our extended valence meson wave function. In Sec. 
Ill we analyze the quark distribution amplitudes and structure functions for the low lying 
mesons and present numerical results. Finally, we discuss and summarize our major results 
in Sec. IV 



II. VALENCE QUARK MESON WAVE FUNCTION 

Following our previous paper we specify the quantization surface S by a timelike four 
vector with n 2 = 1, although the analysis is also appropriate for the case when is 
null-like, n 2 = 0. The transverse and longitudinal components of an arbitrary four vector 
A^ are denoted as A T and A L , respectively, A^ = (A L , A T ) and are defined by 

A L = n ■ A , A T = N 1 - A L n^. 

Since n ■ At = 0, A T has only three independent components which we label (A^, A T ) and 
henceforth identify by At- We define the wave function ^> a {kTu \] Tj), i = 1 ... N as the 
probability amplitude for finding N constituents (quarks, antiquarks, gluons) with transverse 
momenta kxi, helicities A, and flavor-color components Tj in a meson state a, (a = 1 ... 8 for 
the pseudoscalar octet) with momentum Pt = J2i kn, Pl = \j 1 M 2 + \P%\ by the following 
matrix element 

(k Tl ; K; nlPr, a) = {2irf ^2P L 5\P T - £ k Tl ) 

i 

Here the single particle states \kx\ A; r) describe effective, massive constituents quantized on 




(2.1) 
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the spacelike surface £ perpendicular to n M and the longitudinal momenta are constrained 
by the on shell condition, 



kn = \/mf + \kf T \. 



The single particle states are normalized according to 

(k' T ,\',T , \k T ,\,T) = ^{2ir) s 6 3 {k' T -k T )6 x , x 6 T , T 
m 

for fermions and 



(2.2) 



(2.3) 



(P T ,a\Fr,P) = 2{2nfP L 8\P! r - P T )6, 



a/3 



(2.4) 



for bosons. 

The set of wave functions defined by Eq. (|2.1|) constitutes a Lorentz group representation 
basis. In a formalism with a fixed quantization surface, Lorentz transformations depending 
on interactions do not leave £ invariant interactions do not conserve particle number and 



therefore as the system evolves different Fock sectors mix fL4| . This in general leads to 
complicated transformation properties for the wave functions under the action of interaction 
dependent generators of the Lorentz group. A simplification is usually made by using an 
interaction free transformation rule for the wave function in a given Fock sector ||15|| . The 
valence qq wave functions describing a meson state with 4-momenta P and P' = AP = 
C(P — > AP)P respectively are thus related by 



n 




*°(fcri;A i ;r < ) = IIE £) A i v.(i2w) 



1 



i A 



2P' L 




k' 

U\l- 



^{k'^Kin) (2.5) 



with k[ = Aki and Dxy{Rw) denoting the matrix representations of the Wigner rotation 
Rw in the spin space expressed in terms of qq momentum variables, 



R 



R w (k iT ) = C{AP -> P)C{P -> AP)C{P -> P), 



and P = (M, 0) being the meson rest frame momentum. Using this approximation the meson 
state is defined in a standard way and contains only qq valence component. Eqs. fl2.1p, ( |2.8| ) 
and then lead to the following wave function normalization 



]T J [dk lT ]^ a (k tT ; A; r)^(k lT ; A; r) = 5 Q/3 , 



[dki T ]p = [dk iT } P 



2 rl'^b 2 

' A ' T (27r) 3 <5 3 (P T -$> T ). 



n 



(2vr) : 



(2.6) 



In general, approximations and specifically Fock space truncations destroy covariance. Here 
the truncation generates noncovariance by the emergence of an unphysical dependence of 
matrix elements upon n M . In M we describe a method for restoring covariance by allowing 
the quantization surface S, or equivalently the quantization vector 77/ 1 , to transform actively 
under Lorentz transformations by relating n fl to the meson external momenta. 

The normalization Eq. ( j2.6|) is identical to the corresponding nonrelativistic expression 
in the meson rest frame, since the relativistic wave function is constructed to reduce to the 
nonrelativistic one in this frame. For the ground sate pseudoscalar octet the form of the 



relativistic wave function is then derived from Eq. (2.5) to be 



At 



-T1T2 

A a / 
V2 V3 



i{k iT \ Xi) = V2 v =u(k 1T ; Ai)7 5 w(fc 2 T; A 2 ) 

' M. 2 — {mi — m-i) 2 



(2.7) 

where A a are the Gell-Mann SU(3) flavor matrices, I is the identity matrix in the color 
space, mx, m.2 are the quark and antiquark constituent masses respectively, M. is the qq 
invariant mass, 



M{h) = {h + k 2 f = (MM + MM) 2 + {kiT + k 2T ) 



(2.8) 

with kiiikx) given be Eq. ( [2.2| ) and $(A^) being the spin independent orbital wave function 
usually assumed gaussian, 



$(M) = Nexp 



(2.9) 



The overall normalization constant N is determined from Eq. ( p.6| ). 

As explained above the explicit form of the transverse variables will depend on the choice 
for n M which in turn will be specified after selecting which matrix element is to be calculated 



with the wave function of Eq. ( |2.7|) . The wave function specified by Eq. (|2.7|) has also 
been extensively studied in [|] for fixed quantization schemes. Here we study the extension 
of Eq. (|2.7| ) to the most general Dirac structure for the qq system. We write the spinor 
component, £(kiT] Aj) of the wave function in the general form, 

i{k iT - A) = 5>(fcri, Ai)I>(fc T2 , A 2 ). (2.10) 
v 

with the Lorentz 4x4 matrices T p represented by combinations of Dirac matrices and the 
constituent momenta, fcj. Using the Dirac equations for the free u and v spinors it is easily 
shown that sum in Eq. fl2.10|) reduces to two terms involving either 75 or [$i,$2]75- The 



most general wave function \I/ Q for a pseudoscalar meson can thus be written in the form 



v/-M 2 — (rrii ~ rn>2) 2 
u{k TU Ai) [l^p{M) + $1, Ui^t{M)} v(k T2 , A 2 ). (2.11) 



We shall assume that the two spin independent wave functions appearing in Eq. ( J2.ll ) 



can be parameterized with a gaussian shape of Eq. ( |2.9| ) having the same momentum size 
parameter f3 

$ P (M) = ${M) , §t{M) = r PT j^$(M). (2.12) 

Since the tensor term in Eq. (|2.11|) explicitly involves higher powers of the constituent 
longitudinal and transverse variables whose ranges are related to M. and /3, respectively, we 
have chosen the Aif3 factor as the relative normalization between the $p and $t terms in 
Eq. ( |2.12| ). The dimensionless, numerical coeficient rpx will then be determined by fitting 
various properties of the pion. In the rpp = limit the wave function of Eq. (|2.11| ) reduces 
to the one of Eq. ( |2.7| ) and corresponds to the S = L = state in the meson rest frame. The 
tensor term introduces an S = L = 1 orbital component which mixes the S = L = 1 lower 
with S = L = upper components of the Dirac spinors to give an overall J p = 0~ state. 



III. NUMERICAL RESULTS 
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A. 7T electromagnetic form factor 



Before detailing the quark distributions given by the generalized wave functions including 
the additional tensor component $t we shall present new results for the pion electromagnetic 
form factor. The form factor calculation for both noncovariant and covariant variable front 
models has been previously summarized H]. For simplicity we maintain the same quark 
masses and oscillator size parameter (3 used in our previous calculations with m q = (3 = 
250 MeV and vary rp? to optimize the form factor description. It is significant to note 
that the tensor term drastically modifies the form factor behaviour especially in the high 
momentum transfer region Q 2 > 2 GeV 2 . Further as shown in Fig. 1 the generalized model 
provides an agreement with the data in this region without altering the correct low Q 2 
behaviour. As shown below the importance of the tensor term is further demonstrated in 
the analysis of the distribution amplitude and the structure functions. 



B. Distribution amplitude 

The quark distribution amplitude for a meson a, <f) a {!;), can in principle be obtained 
from the moments (£ n ), H 

(O = / dtCFiO (3.1) 

which are formally defined by the matrix elements 

{0|^(0)^75 y ■ ■ ■ d, n ^(0)|P; a) = ifaiO^pP^ ■ ■ ■ ^ + trace (3.2) 

where the trace terms corresponding to higher twist operators have been omited. Here f a 
is the meson decay constant and the normalization is such that = 1. The expansion of 
the quark field operators, ip, ip in terms of constituent qq creation and annihilation operators 
Eqs. ( EZXq ), ( fOg ) leads to 



VQ f d 3 q T 1 m x m 2 / 2g| k 1L - k 2L 

f a J (2tt) 3 JW a Jk 1L k 2L M 2 - (mi - m 2 ) 2 ) V M « M « , 



k 1L m 2 + k 2L m 1 2((k 1L + k 2L ) 2 - (mi + m 2 ) 2 ) ( (mi - m 2 )(k 1L m 2 - k 2L mi)\ 

r PT 1 

m 1 m 2 ' /3(mi + m 2 ) \ (k 1L + k 2L )m 1 m 2 J 



(3.3) 



with q T = (g3,g^), 



k iL = ymf + \q%\, 

M = k 1L + k 2L , (3.4) 

and M a being the meson mass. The decay constant f a is calculated from Eq. ( |3.3| ) using 
= 1 (see also Sec.V in Ref. (§]). In QCD, it can be shown that the physical part of the 
distribution amplitude a (£) is restricted for — 1 < £ < 1 || and for large n the moments 
in Eq. ( |3.2|) behave as (£ n ) — > 1/n 2 implying that a (£ — > ±1) — > 0. In order to reproduce 
this feature in our model calculation we make the replacement 



M a — > M. = kn + A; 2L = ^/m 2 + |g 2 | + ^m 2 + |g 2 |. (3.5) 

In Table 1 we list values for both 7r and K mesons decay constants calculated from Eq. ( |3.3|) 
using spin averaged meson masses M x = 610 MeV and Mk = 790 MeV and compare with 
results using dynamical masses for M a determined by Eq. ( |3.5|) . Note that the experimental 
values lie almost midway between the two methods. The sensitivity to the mass prescrip- 
tion for normalized quantities like moments of the distribution amplitude or structure 
functions, which we analyze in the following subsection, is even smaller. Using Eq. ( |3.5| ) in 
Eq. ( |3.3|) we make the following change of variables 

3 t( \ — 2q T , kiL - k 2L 
"> e(«r) = m + — Jj— (3-6) 

For fixed q^ the variable £,(qr) maps the domain of the variable q^, — oo < gf. < oo into the 

finite interval — 1 < £(qr) < 1- The expression for the distribution amplitude can now be 

obtained from Eqs. (|3TT|), (|3.3| ), ( |3.5|) and ( |3.6|) and is given by 



fa J (2tt) 3 Jk 1L k 2L M 2 - (mi - m 2 ) 2 ) 



ki L m 2 + k 2L mi 2((ki L + k 2L f - (mi + m 2 ) 2 ) ( (mi - m 2 ){ki L m 2 - k 2L mi)\ 

m 1 m 2 /3(mi + m 2 ) \ (k 1L + k 2L )m 1 m 2 J 

(3.7) 
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with kiL and A4, now functions of obtained from Eq. Q3.4Q using Eq. ( [3.6|) to express 
<?t = 9t(0- J{QtiO i s the Jacobian of the transformation from to (£,5r)- The 

explicit forms of the functions, kn^q^,^), ,£) and J(?x>0 are given in the Appendix. 

In Figs. 2 and 3 we plot the function for 7r and mesons respectively. In Fig. 2 
the solid line gives our prediction for the set of parameters m q = (3 = 250 MeV, rp? = 1.0 
which best describe the form factor shown in Fig. 1. The dashed line is the prediction for 
the pion distribution amplitude without the tensor term in the wave function of Eq. ( [2.11D , 
i.e. with rpx = 0. The wide, camel like shape of the distribution amplitude obtained for 
rpT = 1.0, a value that optimizes the form factor description, is similar to that provided by 
the QCD sum rule approach || for the matrix element of Eq. ( |3.2|) , however, the predictions 
for the lowest moments (£ 2 ) = 0.27, (£ 4 ) = 0.13, are smaller then the ones of Chernyak and 
Zhitnitsky ((£ 2 )cz = 0.40, (£ 4 )cz = 0.24) yet closer to those obtained in Ref. |7j and lattice 
calculations [|,f|. The asymmetry in the K distribution amplitude in Fig. 3 is due to the 
large SU(3) breaking due strange quark mass (m s = 480 MeV). Notice the sensitivity of 
the pion distribution amplitude to Tp? which is shown in Fig. 4. The largest sensitivity is 
observed for Tp? ~ 0.5 where the interference between the pseudoscalar and tensor terms 
dominate. 



C. Structure functions 

Structure functions contain important hadronic information and are obtained from in- 
elastic processes. Accordingly we wish to further test our approach by computing the meson 
structure functions and comparing with available data usually extracted from Drell-Yan lep- 
ton ||T0"| , |iTf or charged hadron production |l^] on nuclear targets using mesonic beams. The 



extraction from hadron production experiments is typically more complicated and somewhat 
model dependent due to uncertainties in the hadronization mechanism. In this paper the 
experimental pion structure function was extracted from the Drell-Yan muons produced by 
252 GeV pions on tungsten JTT]]. The theoretical cross section for muon production as a 
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function of longitudinal momentum fraction £p of the muon pair is given by 

d 2 a R ^a 2 \fi( Xl )G»(x 2 ) + tf(*i)G? (x a )] 



dg 2 de F 9Q 4 (eF + 4Q 2 /^) 1/2 

x li2 = [±f F + + 4Q 2 / S ) 1 / 2 ]/2, (3.8) 

where Q 2 an d s are t ne mass of the muon pair and the square of the cm energy respectively, 
and a is the electromagnetic fine structure constant. fZ s )(%i) is the pion valence (sea) quark 
structure function and G^(x 2 ) parameterizes the nuclear contribution. Assuming isospin 
symmetry for 7r~ we have 

f£(x) = xu(x) = xd(x), 

fs(x) = xu s (x) = xu s (x) = . . . = xs(x). (3.9) 

The nuclear contributions can be similarly parameterized in terms of valence and sea quark 
distributions of individual nucleons. The normalization (K factor) is measured to be K ~ 
1.75 ± 0.13 while a perturbative analysis to first order in a s gives K ~ 1.4. The structure 
functions are in principle functions of of Xi and Q 2 , however, since we are describing the 
average data for 36.0 < Q 2 < 72.3 GeV 2 we have suppressed the explicit Q 2 dependence. 
The meson structure functions can equivalently be defined in terms of a diagonal matrix 



element involving the commutator of two vector currents [[16 



Wt u (x,Q 2 ) = i- I dze^{p\^ a ^)^il v ^m\\p) (3-10) 



Arc 

with x = p-q/M 2 , Q 2 = —q 2 > and i referring to a particular flavor. Since it is known that 
the total longitudinal momentum of the hadron is only partially distributed among quarks, 
the constituent quarks which are assumed to cary the entire momentum of the hadron 
cannot be identified with the partons contributing to deep inelastic structure functions. In 
the scaling limit, Q 2 — > oo, Eq. (|3.10| ) can be calculated using a short distance expansion of 
the bilocal operator. In perturbative QCD scaling violations can be described in terms of 
a convolution of the partonic distributions defined at a scale of reference Qq < oo and the 
Altarelli-Parisi splitting functions which characterize the single parton response amplitudes 
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for the change of scale due to radiation of gluons p7|] . Here we also use the convolution 
approach |18|] to relate our constituent quark model to the parton model. In a convolution 
model the quark distribution function for a hadron a, gf , is represented by a product of 
the distribution function of a constituent quark, Q", in a hadron and the probability for a 
constituent quark to fragment into a QCD parton i, q^, 

q«(x,Q 2 ) = J^Q?(x,Q 2 )q i/v (-,Q 2 /Ql). (3.11) 
The constituent quark distributions are defined through Eq. ( |3.10|) with the QCD fields 



replaced by an effective constituent quark/gluon basis at ~ 1 GeV 2 . The Q 2 evolution 
of qi/ v (x/y,Q 2 /Ql) is governed by perturbative QCD. However for any value of Q 2 phe- 
nomenological input is still required. For the average Q 2 ~ 50 GeV of the Drell-Yan data 
the valence and sea quark distributions are 

J* y y 

f 2 -Q n Mis /v (-). (3.12) 
J* y y 

For the valence quark contributions in n~ , q v = q d = q^, while for the sea q s = qd = qj = 
■ • ■ q$. The number and momentum sum rules are respectively 

i 

dxql(x) = 1, 



o 



2 f 1 dxxql(x) +6 f 1 dxxq^(x) = l-g v , (3.13) 
Jo Jo 

where g n represents the fraction of the gluon momentum in the pion currently measured as 



~ 0.47 ||rTJ. The parton distributions in a constituent quark, are usually normalized 
to describe the relevant features of the low-x hadron scattering phenomenology. The Regge 
behaviour at small x motivates the following parameterization 



r 




1 


r| 




T(A) V%~ 



g s/ ,(x) = -(l-x) D - 1 (3.14) 

x 
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with the parameters A, D and C constrained by the momentum sum rule of Eq. (|3.13|) , 

1 C 

+ 6— = 1 - g v 



2A + 1 D 

Each can be determined by comparing the calculated quark distributions and to the 
data. The matrix element which determines the constituent quark distribution functions 
obtained from the light cone expansion of the current product in Eq. (|3.10|) has a very 
similar structure to Eq. ( j3T2|) that defines the quark distribution amplitude. Using the 
techniques developed in the previous section the following expression for Q„(x) is easily 
derived, 



2„_L .,■■>, >91 2 



M 2 - (mi + m 2 ^ 2 

1 + r PT -TT-. ■ 

pM 



(3.15) 



For a = tt~ , mi = m 2 the integrand is a symmetric function of £ = 2x — 1 yielding 
Qjl ( x ) = Qd ( x )- For the kaon, m 2 = m s > mi = m u = m d and the nonstrange and 
strange quark distributions correspond to £ = 2x — 1 and £ = —2x + 1 respectively. Again 
M. = M.{cLt-,£,) an d J(Qt,0 are specified in the Appendix. 

In Fig. 5 we plot valence quark structure functions for the pion calculated in our model 
with tpt = 1.0 (solid line) obtained from the form factor fit and the result without the tensor 
component having Tp? = (dashed line) and compare with data. The comparison suggests 
the importance of the tensor term at large x, 0.6 < x < 0.9 where the sensitivity to the 
constituent quark distributions is the highest. In Fig. 6 we also compare our results for the 
sea quark distributions (solid line) with the curve used for experimental fitting (dashed line). 
The value C = 0.086 for the parameter in Eq. ( ft.14 ) is obtained by requiring the theoretical 



and experimental curves to agree at x — 0. Fitting A by the valence quark distribution 
yields A = 0.75 and the sum rule of Eq. ( |3.13| ) then gives D = 4.0. These numbers are in a 
good agreement with the ones obtained from the unpolarized nucleon structure function fits 



confirming hadron independence of the splitting functions qi/ v (x) [|T^,|T^]. In Fig. 7 we also 
show our predictions for the strange (solid line) and light quark structure functions (dashed 
line) in kaon. 
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IV. SUMMARY AND CONCLUSIONS 



Within the framework of our covariant variable front quark model, we have generalized 
the constituent, qq pion wave function and have studied the distribution amplitudes and 
structure functions. The extended model provides excellent agreement with the experimen- 
tal data for the structure functions and the electromagnetic form factor as well as a good 
description of the decay constant. The improved description is due to a tensor component 
in the wave function which is a relativistic correction in the rest frame. We have computed 
the structure functions and the distribution amplitudes with the wave function of Eq. ( |2.11fl 
in a noncovariant light cone quantization scheme and find results are similar to the vari- 
able front model. This confirms front independence of our results and is consistent with a 
previous assertion 0] that different front formulations do not lead to significant differences 
in the predictions for various mesonic properties. The magnitude of the tensor term which 
optimizes both the form factor and structure function now indicate a camel shape for the 
distribution amplitude with the dip at £ = although not as profound as the one suggested 
by the old QCD sum rules, but quite similar to recent nonlocal QCD sum rule calculations 
and lattice results. Since the shape and the moments of the quark distribution amplitude 
are very sensitive to the interference between the pseudoscalar and tensor components of the 
wave function the detailed knowledge of former will provide important information on the 
structure of the meson wave function. As mentioned in the introduction, the distribution 
amplitude is not directly measurable, however, information on the magnitude of the lowest 
moments can be extracted from the meson form factor describing the hadronic part of the 
e meson — > e~y transitions [§,0]- Finally, distribution amplitude studies provide an effec- 
tive forum for theoretical comparisons of alternative model formulations which can provide 
significant insight into QCD dynamics. 
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V. APPENDIX 



The variable change of Eq. ( |3.6| ) gives, 

M?(^)(l-0 2 -M 2 2 (^)(l + 2 



/8(i-a^i(^)(i-o+^(^)(i+o] 

where 



£{<1t) = < + (?t) 



^(^,0 = V^ 2 (9t) + (#(^0) 2 > 



and the Jacobian J(qx,0 * s given by 



J(#, = 2— 2-5^ + (A^ - 2^Afj] - -AfMjI - 2£ 



1-e 2 l-^ 2 s /I 2 + (A/i) 2 -2e7lA y u 



+2( i -e) 7z3(A/i)3 

with 



A/i = A//(g£) = -{hMt) ~ friOr))- 
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FIGURES 



FIG. la,b. Pion electromagnetic form factor. Solid line shows the result obtained with 
the qq wave function of Eq. (2.12) for r PT = 1.0. Dashed line is the result of the covariant 
model of Ref. [4], dash-doted and doted lines are the results of relativistic models from 
Refs. [2] and [3] respectively. Data is taken from Ref. [21]. 

FIG. 2. Pion distribution amplitude with {r-pp = 1, solid curve) and without {r-pp = 0, 
dashed curve) the contribution from the tensor term in the qq wave function. 

FIG. 3. Same as Fig. 2 for the K distribution amplitude. 

FIG. 4. rpT dependence of the pion distribution amplitude. 

FIG. 5. Valence quark pion structure function f£(x) = xq v (x) (curve convention same 
as in FIG. 2.). Data is taken from Ref. [11]. 

FIG. 6. Sea quark pion structure function. Solid line gives our model prediction with 
r PT = 1., dashed line shows experimental parameterization, f£(x) = 0.173(1 — a;) 8 - 4 of 
Ref. [11]. 

FIG. 7. Valence light quark (dashed line) and strange quark (solid line) kaon structure 
functions. 
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TABLES 

TABLE I. 7r and K mesons decay constant calculated with the formula of Eq. (3.3). The two 
set of results correspond to the use of spin averaged, constituent meson masses M and dynamical 
masses M respectively, (m u = m d = 250MeV, m s = 480MeV, (3 = 250MeV, r PT = 1.0) 



M 


U[MeV] 
101. 


f K [MeV] 
137. 


M 


75. 


104. 


exp. 


93.2 


113. 
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